The long activations of a2 glycine channels can be described 
by a mechanism with reaction intermediates ("flip") 
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The a2 glycine receptor (GlyR) subunit, abundant in embryonic neurons, is replaced by al in the adult nervous 
system. The single-channel activity of homomeric ot2 channels differs from that of al-containing GlyRs, as even at 
the lowest glycine concentration (20 uM), openings occurred in long (>300-ms) groups with high open probability 
(fopen! 0.96; cell-attached recordings, HEK-expressed channels) . Shut-time intervals within groups of openings 
were dominated by short shuttings of 5-10 us. The lack of concentration dependence in the groups of openings 
suggests that they represent single activations, separated by very long shut times at low concentrations. Several pu- 
tative mechanisms were fitted by maximizing the likelihood of the entire sequence of open and shut times, with 
exact missed-events allowance (program hjcfit) . Records obtained at several glycine concentrations were fitted 
simultaneously. The adequacy of the different schemes was judged by the accuracy with which they predicted not 
only single-channel data but also the time course and concentration dependence of macroscopic responses elic- 
ited by rapid glycine applications to outside-out patches. The data were adequately described only with schemes in- 
corporating a reaction intermediate in the activation, and the best was a flip mechanism with two binding sites and 
one open state. Fits with this mechanism showed that for a2 channels, the opening rate constant is very fast, 
^130,000 s -1 , much as for alp GlyRs (the receptor in mature synapses), but the estimated true mean open time is 
20 times longer (around 3 ms) . The efficacy for the flipping step and the binding affinity were lower for a2 than 
for al(3 channels, but the overall efficacies were similar. As we previously showed for al homomeric receptors, in 
a2 glycine channels, maximum P open is achieved when fewer than all five of the putative binding sites in the pen- 
tamer are occupied by glycine. 



INTRODUCTION 

Out of the four known isoforms of the glycine receptor 
(GlyR) a subunit, a2 (Grenningloh et al., 1990; Akagi 
et al., 1991) may play a major role in the embryonic 
form of GlyRs (Legendre, 2001), as its transcripts, abun- 
dant in the embryonic and neonatal rat, are gradually 
replaced by al in the adult (Akagi et al., 1991; Malosio 
et al., 1991; Watanabe and Akagi, 1995). This is sup- 
ported by a variety of functional observations. Takahashi 
et al. (1992) observed longer apparent open times in 
recombinant a2 versus al GlyRs, and in embryonic 
versus adult rat spinal GlyRs. The decay of glycinergic 
IPSCs also speeds up with development. The simplest 
explanation is that embryonic and neonatal GlyRs con- 
tain mostly a2 subunits, whereas adult receptors are 
al(3 heteromers (Becker et al., 1988; Langosch et al., 
1988; Hoch et al., 1989) . Nevertheless, it is far from cer- 
tain that the slow glycinergic synaptic currents observed 
early in development are mediated by a2 homomeric 
receptors. Homomeric GlyRs probably cannot cluster at 
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the synapse, because only (3 subunits can interact with 
gephyrin (Meyer et al., 1995). Second, a2 is expressed 
in areas where no glycine-mediated synaptic transmis- 
sion has been detected (Flint et al., 1998; Mangin et al., 
2002) and by nonneuronal cells, such as glial or pro- 
genitor cells (Belachewetal., 1998; Nguyen et al., 2002) . 
Mangin et al. (2003) showed that a2 homomeric recep- 
tors are unsuitable for fast synaptic transmission because 
they activate and deactivate slowly and have a low prob- 
ability of opening in response to the rapid application 
of brief glycine pulses (Mangin et al., 2003) . It has there- 
fore been suggested that a2 homomeric receptors may 
be mainly extrasynaptic and that they can be activated 
by a paracrine-like release of glycine or taurine during 
neuronal development (Legendre, 2001; Clements, 2002; 
Mangin et al., 2003). In addition to that, synaptic cur- 
rents recorded in neonatal hypoglossal motoneurons 
decay almost 10 times faster than currents elicited on 
recombinant a2 channels (Mangin et al., 2003), lead- 
ing Singer et al. (1998) to suggest that they are medi- 
ated by a23 heteromeric GlyR. 
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Global fitting of kinetic models to single-channel data 
from recombinant al-containing GlyRs led us to de- 
scribe an activation mechanism ("flip") that explicitly 
incorporates intermediate shut states between agonist 
binding and channel opening (Burzomato et al., 2004) . 
Mechanisms that include pre-open shut states such as 
flip states, or for muscle nicotinic receptors, "primed" 
states (Mukhtasimova et al., 2009), are attractive because 
they offer a suggestive kinetic parallel to the events that 
are likely to happen in the extracellular domain before 
channel opening, such as domain closing and/or Gloop 
capping. In addition to that, they have been found to 
describe well channel activation for a variety of wild- 
type and mutant subunit combinations for glycine, 
muscle nicotinic, and GABA A receptors, notably with 
respect to the differences between full and partial ago- 
nists in the nicotinic superfamily and to the complex 
structure of unliganded openings in gain-of-function 
mutants (Burzomato et al., 2004; Plested et al., 2007; 
Lape et al., 2008; Mukhtasimova et al., 2009; Keramidas 
and Harrison, 2010). Our aim in the present work was 
to extend this approach to recombinant a2 GlyRs, to 
establish a quantitative activation mechanism for this 
channel, and to test whether a flip-type mechanism 
could provide an adequate description for both single- 
channel and macroscopic data from recombinant a2 
homomers. For the first time in fitting GlyR channel ac- 
tivity, we found that a flip-type mechanism describes the 
data somewhat better than a less constrained mecha- 
nism incorporating distal shut states (such as the mecha- 
nism proposed by Jones and Westbrook, 1995). 

MATERIALS AND METHODS 

GlyR a2 subunit cDNA construct, HEK293 cell culture, 
and transfection 

HEK293 cells (American Type Culture Collection and LGC Stan- 
dards) were maintained in a humidified incubator at 37°C (95% 
air/5% CO2) in Dulbecco's modified Eagle's medium supple- 
mented with: sodium pyruvate (0.11 g/liter), heat-inactivated 
fetal bovine serum (10% vol/vol), and penicillin G (100 U/ml)/ 
streptomycin sulfate (100 ug/ml; all from Invitrogen). Cells were 
passaged every 2-3 d, up to 30 times. Cells were plated on glass 
coverslips (13-mm diameter, sterilized over flame), placed in 
35-mm sterile Petri dishes containing 1.5-2 ml Dulbecco's modi- 
fied Eagle's medium, 4—6 h before transfection, to allow cells to 
attach. Cells used for fast concentration jumps were plated onto 
coverslips coated with poly-L-lysine (Sigma-Aldrich) by incuba- 
tion of the coverslips for ~40 min, followed by washing with auto- 
claved water and oven drying at 60°C. 

HEK293 cells were transiently transfected with the rat GlyR a2 
cDNA (Isoform A; available from GenBank/EMBL/DDBJ under 
accession no. AJ310837; full-length sequence verified) in the 
pcDNA3.1 vector. Transfection was by calcium phosphate-DNA 
coprecipitation (Groot-Kormelink et al., 2002). In brief, cDNA 
was added to CaCl 2 solution (340 mM in sterile water) at a volume 
ratio of 1:5, and this mixture was added drop-by-drop to an equal 
volume of 2* Hank's buffered saline (280 mM NaCl, 2.8 mM 
Na 2 HP0 4) and 50 mM HEPES, pH 7.2 with NaOH) to form the 
precipitate. The final mixture was then added dropwise over the 



plated cells. The cDNA used for transfection contained 3-20% 
a2/pcDNA3.1, 18% pEGFP-cl, and 62-79% pcDNA3.1 vector 
lacking the a2 coding sequence. Adding empty vector was found 
to be important for the health of the transfected cells and to 
ensure that the level of a2 GlyR expression was optimal for re- 
cording. The total amount of the final cDNA mixture per plate 
was kept constant at 3 ug. Recordings were performed 2-3 h after 
washing and no more than .3 d after transfection. We attempted to 
obtain single-channel recordings also from a2|3 heteromeric 
GlyR, but the expression level was consistently too low, despite 
considerable effort to improve it. 

Single-channel recording 

Cell-attached single-channel currents were recorded at 19-21°C 
with thick-walled borosilicate pipettes (with filament; Harvard Ap- 
paratus) coated near the tip with Sylgard (Dow Corning) and fire- 
polished before use to a final resistance of 5-12 MO. Electrodes 
were filled with extracellular medium (in mM: 20 Na-gluconate, 
102.7 NaCl, 2 KC1, 2 CaCl 2 , 1.2 MgCl 2 , 10 HEPES, 20 TEA-C1, 
15 sucrose, and 14 glucose, pH 7.4 with NaOH; osmolarity 320 
mOsm) plus the appropriate concentration of glycine. All solu- 
tions were prepared in HPLC-grade water (VWR International) to 
minimize contamination by glycine and filtered through a 0.2-pm 
Cyclopore track-etched membrane (GE Healthcare) to remove 
impurities that can block electrodes or affect the quality of the 
seal. The bath level was kept low to reduce noise. 

Currents were recorded with an amplifier (Axopatch 200B; 
MDS Analytical Technologies) at a pipette potential of +100 mV, 
with no correction for junction potential, as this is calculated 
to be at most +1 mV (Clampex 9; MDS Analytical Technologies). 
Single-channel currents were prefiltered at 10 kHz by the amplifi- 
er's four-pole Bessel filter and saved on a digital audio tape with a 
recorder (DRA-200; Bio-Logic). Recordings were replayed, filtered 
at 5 or 7 kHz, and sampled at 50 or 71.43 kHz, respectively, depend- 
ing on the signal-to-noise ratio, with a Digidata 1 322A and Clampex 
software (MDS Analytical Technologies). 

Single-channel current analysis 

Single-channel currents were idealized by time-course fitting 
(SCAN program). Segments showing drifts in the baseline, seal 
breakdowns, or dubious channel openings were excluded, and 
adjacent intervals were marked unusable as appropriate. After 
imposing a resolution of 20-25 ps (for both shut and open times) , 
the number of resolved time intervals (sum of shut-time and 
open-period intervals) ranged from 498 to 2,368 in the 12 patches 
used in the maximum likelihood fits. 

Idealized current amplitudes longer than twice the filter rise 
time (66 ps at 5 kHz) were plotted as an amplitude histogram 
(EKDIST program). The main conductance level accounted for 
>99% of all fitted open intervals, so we pooled together all open- 
ings for kinetic analysis. The average amplitude of openings to 
the main conductance level was 5.9 ± 0.2 pA (n = 25 patches) but 
varied considerably, ranging from 4.5 to 9.1 pA. In the 12 patches 
used in our analysis, the range of mean amplitudes was 4.9-8.0 pA 
(6.1 ± 0.3 pA). This variability in amplitude was expected for cell- 
attached records because it will be affected by the intracellular 
chloride concentration and the resting membrane potential, both 
of which vary from cell to cell. Variability in the intracellular chlo- 
ride concentration and in the true transmembrane potential is a 
source of error in so far as it affects the rate constant values in re- 
ceptor activation. 

The distributions of apparent shut times and open periods 
were fitted with a mixture of exponential probability density func- 
tions. Each open period was defined as the duration of time for 
which the channel remained continuously open, regardless of 
amplitude; that is, the time between two adjacent shut times that 
were both longer than the imposed resolution. 
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The true number of the channels in each patch is unknown, but 
for kinetic analysis, it is essential to isolate stretches of openings 
that are likely to arise from one individual channel. This was done 
in the initial stage of the analysis by fitting shut-time distributions 
with EKDIST to determine a critical shut time (t clk ) that would 
allow us to divide the record into true one-channel stretches of 
openings. Fits to the distributions obtained at EKDIST are only de- 
scriptive and not mechanism dependent; hence, parameters esti- 
mated by EKDIST were not used further in the analysis. The criterion 
for a group of openings to be judged to have come from one indi- 
vidual channel is essentially that the probability of being open 
within the groups is high enough that we can be almost certain that 
a double opening would be seen if more than one channel were 
present (Colquhoun and Hawkes, 1990). Groups with fewer than 
10 open periods were excluded. For fitting, we tested different val- 
ues of trit (with or without the use of Colquhoun-Hawkes-Srodzinski 
vectors; see below) . We found that the most consistent fits were ob- 
tained when we used the Z crit values reported in the Results. 

For kinetic analysis, the postulated mechanisms were evaluated 
using maximum likelihood fitting and the program HJCFIT. For 
this purpose, the idealized data from all concentrations were 
grouped into three independent sets, with each set containing 
four patches at different agonist concentrations (20 or 30 uM, 
100 |iM, and 1 and 10 mM glycine). All data from one set were 
input into the program to be fitted simultaneously to a specified 
mechanism, together with the resolution of 20-25 ps (the dura- 
tion of the fastest event that can be detected) , a f CI i t value and the 
glycine concentration for each patch. The program calculates a 
likelihood value for the dataset from the initial guesses for the 
rate constant values, taking into account the imposed resolution, 
which is used to correct for missed events (Hawkes et al., 1990, 
1992) . The rate constant values are then adjusted to maximize the 
likelihood until a maximum is reached. 

Shut times longer than t clit cannot be used directly in the fitting 
process because they are affected by the (unknowable) number of 
channels in the patch. Nevertheless, some information is recover- 
able from the long shut times. The true (one-channel) shut time 
between one activation and the next is not known, but it is known 
that it must be at least as long as * cr j t . This information can be used 
to improve the fit by replacing the normal steady-state vectors as 
the start and end of each group of openings, with the Colquhoun- 
Hawkes-Srodzinski vectors described by Colquhoun et al. (1996; 
Eq. 5.11). This was done for low concentration records (in this 
case, all concentrations except 10 mM), and simulations show 
that it is useful (Colquhoun et al., 2003) . It cannot be used at con- 
centrations that are high enough that desensitization occurs, 
because it is model dependent and our models do not include 
desensitized states. 

We tested whether our fits were robust to the choice of initial 
guesses for the rate constant values to exclude the possibility of 
convergence to local maxima in the likelihood surface. If the fit is 
good and the likelihood surface has only one well-defined maxi- 
mum, changing the initial guesses should not affect the estimates 
for the rate constants. This was not always the case, and some- 
times several equivalent solutions (with similar likelihood values) 
were possible, suggesting that the likelihood surface was rough 
with several maxima. In this case, the best solution was chosen on 
the basis of fit quality, physical plausibility of the rate constant es- 
timates, and predictive power of the fitted rate constants. For ex- 
ample, we discarded fits that gave physically impossible association 
rate constant values (i.e., >10 9 M _1 s _1 ) . As a test of fit quality, we used 
several data displays: shut-time distributions, open-period distri- 
butions, and the predicted open probability (P open )-concentration 
curve. We also used the results of the fits to simulate realizations 
of single-channel activity (SCSIM program) and to calculate macro- 
scopic current relaxations to a glycine concentration jump (SCALCS 
program) , and we compared them to our experimental observations. 



All analysis programs are available from http://www.ucl.ac.uk/ 
Pharmacology/ dcpr95.html. 

Macroscopic currents: recording and analysis 

Macroscopic currents evoked in outside-out patches by agonist 
concentration pulses were recorded with thick-walled borosilicate 
pipettes (with filament; Harvard Apparatus), fire-polished to a 
final resistance of 8-15 Mfl when filled with "low [CP] solution" 
(in mM): 121.1 K-gluconate, 1 CaCl 2 , 1 MgCl 2 , 10HEPES, 11 EGTA, 
6 TEA-C1, and 2 MgATP, pH 7.2, adjusted with KOH; osmolarity 
325 mOsm, adjusted with sucrose. Because of the very low CP con- 
centration in this solution, the silver/ silver chloride wire of the 
pipette holder was re-chlorided before each experiment to compen- 
sate for the increased flow of CP ions from the wire to the solution. 
The bath solution was the same as for single-channel recordings. 
All outside-out recordings were performed at a pipette holding 
potential of —100 mV. We did not correct for junction potential, 
and this gave a maximum error of +14.6 mV in low [CP] solution. 
During the 3-5-min intervals between glycine applications, patches 
were held at — 60 mV to prolong the stability of the seal. 

Glycine, dissolved in bath solution, was applied to outside-out 
patches with a theta tube (Hilgenberg GmbH) cut to a final diam- 
eter of ^T50 pm at the tip. The tube was driven by a piezo stepper 
(Burleigh Instruments, Inc.). The exchange time was measured 
by the application of 30% diluted bath solution before the experi- 
ment (to optimize the electrode position) and after the rupture 
of the patch. Macroscopic currents were recorded with an Axo- 
patch 200B, digitized with Digidata 1322A, and saved directly on 
a computer via Clampex software (5 kHz filtering, 20-50 kHz 
sampling rate; all from MDS Analytical Technologies) . 

Different glycine concentr ations were applied in random order 
for a duration chosen to allow the current response to reach a pla- 
teau or a peak. All currents were corrected for rundown using a 
standard glycine concentration (about £C 30 ) , which was applied 
every third response. Concentration-response curves were plotted 
using these corrected responses (averaging up to five responses), 
and each experiment was fitted to the Hill equation by least squares 
with equal weights, using the CVFIT program to get estimates of 
maximum current, Hill coefficient ( mh) , and EQ^. Numbers in the 
text are means ± SD of the mean (SDM) from the individual fits of 
each experiment. For the purpose of display, responses were nor- 
malized to the fitted maximum response of the same experiment 
and pooled to give an average concentration-response curve. 

To study the kinetics of macroscopic currents, 10-50 responses 
were recorded for each glycine concentration at intervals of at 
least 20 s to allow full recovery from desensitization (Gentet and 
Clements, 2002). Responses at each concentration for each patch 
were averaged, excluding failures or responses that contained 
patch breakdowns. Only experiments in which the rundown 
between the first and last three responses was <5% were included 
in the analysis. The rise time for average currents and tip poten- 
tials was measured as the time from 20 to 80% of the peak re- 
sponse using Clampfit 9.2 software; experiments in which the 
open tip response had a 20-80% exchange time slower than 
300 ps were rejected. Both the time course of desensitization at 
prolonged glycine applications (200 ms) and the time course of 
deactivation at short applications (2 ms) were fitted with a sum of 
one or more exponentials. 

RESULTS 

Some definitions 

For the purpose of clarity, we shall first define four terms 
that are used throughout to describe patterns of chan- 
nel openings: group, activation, burst, and cluster. 
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The term "group" of openings is used as a generic de- 
scription of any group of openings that is separated by 
short shut times, regardless of why they occur. 

An "activation" is what happens between the start of 
the first opening of a channel until the end of the last 
opening before the receptor returns to the unoccupied 
resting state. The activation is the elementary event pro- 
duced by agonist binding, consists of a burst of one or 
more openings that occur in quick succession, and is 
what underlies the deactivation rate after a macroscopic 
jump to zero concentration (see the appendix of Wyllie 
etal., 1998). 

The term "burst" refers to the groups of openings 
that occur at very low agonist concentrations. At suffi- 
ciently low agonist concentrations, the observed bursts 
of openings are a close approximation to channel acti- 
vations (Edmonds etal., 1995; see Discussion). The shut 
times between activations contain one or more sojourns 
in the unoccupied resting state, as well as bindings that 
happen not to produce an opening. They may also con- 
tain sojourns in desensitized states. 

A "cluster" of openings, defined as in Sakmann et al. 
(1980), refers to the long groups of openings that are 
observed at high agonist concentrations. A cluster con- 
sists of many activations that occur in quick succession 
(because of the high agonist concentration) , so one ac- 
tivation cannot be distinguished from the next one. 
Clusters are separated by shut times that have one or 
more sojourns in long-lived desensitized states. They 
are seen with most ligand-gated ion channels but are 
not easy to identify in a2 receptor records. 



Properties of macroscopic currents from a2 GlyRs 

To understand the ct2 GlyR channel behavior in non- 
equilibrium conditions, we recorded macroscopic cur- 
rents activated by fast application of glycine. We have 
recendy shown (Pitt et al., 2008) that the gating of cd(3 
GlyRs is modulated by internal CI . Channel deactiva- 
tion after short glycine applications becomes faster in 
outside-out patches when intracellular chloride is lowered 
to concentrations similar to those likely to be present 
in the intact HEK cell (i.e., in cell-attached recordings) . 
As this is likely to occur also for ct2 homomeric recep- 
tors, we recorded the macroscopic jumps using a low 
(10-mM) intracellular chloride concentration to make 
sure that in the different experiments, channels were 
exposed to similar chloride concentrations so that the 
information from different patch configurations could 
be combined in our investigation of the activation 
mechanism. Traces in Fig. 1 A are the outside-out cur- 
rent responses of ct2 homomeric GlyRs to the fast appli- 
cation of 2-ms pulses of 10 mM glycine. The results of 
fitting the time course of current responses to short 
pulses of 0.5-10 mM glycine are shown in Table I. Two 
exponential components were needed to fit the current 
decay, and the slower component was the most promi- 
nent by far. 

We also recorded outside-out a2 GlyR macroscopic 
responses to long (200-ms) concentration jumps at 
saturating glycine concentrations (10 mM; —100 mV; 
low [CP] pipette solution) . A typical current example 
is shown in Fig. 1 B. All patches (n = 6) revealed a com- 
plex desensitization phase that could be best fitted with 
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Figure 1. Macroscopic responses of c«2 GlyR to 
rapid glycine application to outside-out patches. 
(A) The current trace represents the average of 
12 sweeps evoked by a 2-ms concentration jump 
of 10 mM glycine applied to an outside-out patch 
(held at -100 mV) from a HEK293 cell express- 
ing a2 GlyRs. Recordings were performed with 
10 mM CP pipette solution. The time course 
of glycine application, taken from the open tip 
response after seal rupture, is shown above the 
trace. (B) Similarly, the average current obtained 
with a long pulse (200 ms) of 10 mM glycine 
(average of 20 sweeps) . Note the complex desen- 
sitization phase (for this patch the decay was fit- 
ted by three components, with time constants and 
areas of 1.8 ms, 0.5%; 25 ms, 1.5%; and 1,300 ms, 
98%). (C) Traces are average currents elicited 
in one outside-out patch by two to five glycine 
applications using 10 mM CP pipette solution. 
Glycine concentrations were applied in random 
order and for a time sufficient to obtain a peak 
current. (D) Concentration-response curve for 
glycine obtained from pooled data from out- 
side-out currents as those in C (refer to Materials 
and methods for details) . Currents were normal- 
ized to the fitted maximum for each experiment, 
and the pooled data were fitted with the Hill 
equation (n = 6). 



200 



Activation mechanism of a2 GlyRs 



mixtures of three exponential components: 1.9 ± 0.3 ms 
(0.6 + 0.2% fractional area) , 23 ± 3 ms (2.2 + 0.6%) , and 
820 ± 180 ms (97.2 ± 0.7%). This complexity of desensi- 
tization was independent of current amplitude: there 
was no apparent relationship between the peak current 
amplitude of the different patches and the number of 
the components needed for fitting currents. 

Finally, to characterize the concentration-response 
relation for a2 channels, we performed concentration 
jumps at a range of glycine concentrations. The traces 
in Fig. 1 C show examples of averaged currents re- 
corded in the outside-out configuration using a low 
[CP] pipette solution (10 mM chloride; -100 mV). 
Concentration-response curves obtained from six 
patches gave an average EC^ of 770 ±100 pM and an n H 
of 1.99 ±0.25 (Fig. 1 D). 

General features of a2 GlyR single-channel activity 
Fig. 2 A shows representative glycine-activated single- 
channel currents recorded in the cell-attached configu- 
ration (pipette holding potential of +100 mV) from 
HEK293 cells transfected with rat a2 GlyR cDNA. It was 
immediately clear that the single-channel activity of a2 
GlyRs is very different from that of al (Beato et al., 
2004) or al(3 GlyRs (Burzomato etal., 2004). Openings 
of al or a 1(3 GlyRs occur in short bursts at low glycine 
concentrations (Fig. 2 B, bottom trace) . These bursts 
are activations that, as glycine concentration is increased, 
become more closely spaced because the time spent in 
unliganded resting states gets shorter until, in high gly- 
cine, they become grouped into long, high P ope n clusters, 
separated by long shut periods. 

In apparent contrast to what is seen with al or a 1(3 
receptors, a2 activity appeared to be very similar at all 
glycine concentrations (Fig. 2, A and B) . Even at 20 pM, 
the lowest glycine concentration at which single-channel 
currents could be detected, groups of a2 channel open- 
ings were long and separated by long shuttings (usually 
lasting for minutes) . Because of the length of the shut 
times, each recording contained a relatively small num- 
ber of these groups of openings and few long shut times. 



Another consequence of the small number of long shut 
times that could be recorded is that it was impossible to 
measure any changes in their length as glycine concen- 
tration was increased. Moreover, the "experimental" 
success rate was generally low, and only a minority of a2 
patches displayed any channel activity (^1 in 5 patches 
at glycine concentrations of 100 pM or greater and only 
~1 in 10 at 20-50 pM) . We could find no single-channel 
activity in records at 10 pM glycine, despite recording 
from >70 good patches for 5-25 min. 

After idealization of single-channel records by time- 
course fitting, the distributions of apparent open and 
shut times were plotted and fitted with mixtures of ex- 
ponential probability density functions. The distribu- 
tions of apparent open periods are shown in Fig. 2 A 
(middle column). These were fitted with up to three 
exponential components (Table II). The longest com- 
ponent was predominant at all concentrations, with an 
average area of 82-98%, and had a time constant of 
^30 ms at all concentrations other than 20 pM glycine, 
where it appeared to be about threefold shorter. Both 
of the two faster components represent a very small 
number of intervals and were poorly defined, and a 
single-exponential density fitted the distribution quite 
well in most cases. The absolute number of apparent open- 
ings or shuttings per patch was small, despite the length 
of the clusters and the length of the recordings. Appar- 
ent shut-time distributions (Fig. 2 A, right column) were 
fitted with a minimum of three components (Table III) . 
At all glycine concentrations, most shuttings were con- 
tributed by the shortest component (area, 85-99%). 
Given its average time constant of 5-10 ps, the vast ma- 
jority of the events in this component were shorter than 
the experimental resolution of 20-25 ps and thus were 
missed. The number of shut-time intervals that fell into 
the two longest components was very small, and conse- 
quendy, these components were not very well defined. 
The time constants of these two components appeared 
to decrease about threefold when glycine concentra- 
tion increased from 20 pM to 10 mM. The existence 
of more than one shut-time component at saturating 
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TABLE I 

Rise-time and decay-time constants of glycine concentration jump (2 ms) current responses 
Amplitude Rise time Deactivation 



20-80% 



t1 

(fractional area) 



t2 

(fractional area) 



mM 


P A 


ms 


ms 


ms 


10 (13) 


170 ±30 


0.42 ± 0.08 


5.6 ± 0.8 


130 ± 10 








(3 ±3%) 


(97 ± 3%) 


1 (13) 


60 ± 10 


1.1 ±0.1 


5 + 2 


92 ±8 








(1 ±2%) 


(99 ± 2%) 


0.5 (9) 


18 ±4 


1.2 ±0.1 


4± 1 


83 + 6 








(1 ± 1%) 


(99 ± 1%) 


The deactivation phase 


was fitted with two 


exponential curves. Recordings were performed 


with low (10 mM) chloride in 


the pipette at -100 mV. Data 



are given as mean + SDM. 
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concentrations (Fig. 2 A and Table III) is not compat- 
ible with mechanisms that have only one shut state 
when fully liganded. Additional fully bound shut 
states are therefore required, either as intermediates 
during channel opening or otherwise distal to the 
binding steps. 

The records shown in Fig. 2 A show a superficial re- 
semblance to the clusters of openings that are seen with 
other receptors at concentrations of agonist that are 
high enough to produce long desensitized intervals, but 
it seems that the appearance is deceptive. These long 
groups of openings appear to be unusually long indi- 
vidual activations of the channel, rather than clusters of 
many activations. This interpretation is confirmed by 
the analysis below. 

The best way to display and characterize the concen- 
tration dependence of single-channel activity is to plot a 



one-channel concentration-P open curve. This is done by 
identifying stretches of openings that are likely to be 
the expression of the activity of a single-channel mole- 
cule by setting for each idealized record a critical shut- 
time duration 4„t value (refer to Materials and methods). 
We attempted to construct a single-channel concentra- 
tion-Popen curve for a2 receptors in this way, but it 
proved to be impossible. A typical display of the -P op en 
values calculated for groups of openings, pooled and 
averaged, is shown in Fig. 2 C: the P op en was high, well 
above 90% at all agonist concentrations. Similar results 
were obtained at all the values of 4„t tested (1-100 ms) . 
This is because (see Discussion) these groups represent 
channel activations, rather than clusters of several acti- 
vations. Hence, the P op en data found in this way yielded 
no information about the concentration dependence 
of channel activity. 
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Figure 2. Single-channel prop- 
erties of a2 GlyRs in the cell- 
attached configuration. (A) Traces 
on the left are continuous 
stretches (20 s) of single-channel 
currents from oi2 GlyRs ex- 
pressed in HEK293 cells (pipette 
holding V, +100 mV) recorded 
in the presence of 20 pM to 
10 mM glycine (5-7-kHz filter). 
On the right are the correspond- 
ing open-period and shut-time 
distributions, fitted with mixtures 
of exponential probability den- 
sity functions. (B) a.2 GlyR single- 
channel currents in 20 pM and 
10 mM glycine are compared 
with equivalent traces recorded 
from HEK293 cells expressing 
al|3 GlyRs, under similar record- 
ing conditions (from Burzomato 
et al., 2004). Note the long 
bursts for tx2 channels at the low 
glycine concentrations. (C) The 
^open concentration-response 
curve for a2 receptor clusters 
at a range of glycine concentra- 
tions is shown. Each point in 
the plot represents the aver- 
age ^open value for all idealized 
patches at each concentration. 
The Popen for each patch was 
calculated as the ratio between 
the total open time and the total 
length of the clusters. 
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TABLE II 





EjflpVflCdl j\t oj JfllXtHTCS oj 


exponential probability density functions to the apparent oj 


icn-period distributions 


[Glycine] 


t1 


t2 


t3 


Mean open time 




Area (%) 


Area (%) 


Area (%) 








WIS 


ms 


WIS 


20 pM n = 2 


0.017 + 0.009 


0.20 


13.1 ± 1.5 


12.5 ± 1.8 




(15 ±10) 


(5) 


(83 ± 12) 




30 pM n = 2 


0.011 ±0.001 


0.52 ±0.13 


37.4 ± 7.3 


34.4 ± 7.5 




(12 ±2) 


(5±1) 


(82 ± 3) 




50 pMn = 2 


0.051 




29.0 ±2.6 


28.7 ±2.9 




(4) 




(98 + 2) 




100 pMn=3 


0.029 


0.562 ± 0.449 


31.0 ±5.1 


30.1 ±4.8 




(3) 


(4 + 1) 


(97 ±1) 




1 mM n = 3 


0.023 ±0.011 




30.4 + 2.4 


29.9 + 2.2 




(11 ±8) 




(93 ± 6) 




10 mM n = 3 


0.018 ± 0.006 




34.8 ± 4.6 


32.9 ± 6.1 




(4.33 ± 0.03) 




(97 ±1) 





n is the number of idealized patches used for each concentration. The time constant value and fractional area (in brackets) of each component are the 
average + SDM for each concentration. 



Maximum likelihood fit of kinetic mechanisms: 
initial considerations 

Next, we performed maximum likelihood fitting of 
kinetic mechanisms to our single-channel data, using 
HJCFIT. Single-channel records were grouped into sets, 
and each set contained four patches that spanned the 
range of glycine concentrations (refer to Materials and 
methods) . Because (as above) we were not able to assess 
the concentration dependence of the single-channel 
-Popen (e.g., Burzomato et al., 2004; Lape et al., 2008), for 
a2 receptors we had to rely on the macroscopic concen- 
tration-response curve (Fig. 1 D) for the choice of con- 
centrations. Ideally, to provide sufficient information 
on the different steps of the mechanism, each set should 
include records from low, intermediate, and saturating 
concentrations of glycine. 



Initial attempts to fit a range of concentrations up to 
1 mM gave inconsistent results. However, the macro- 
scopic concentration-response curve (Fig. 1 D) sug- 
gests that 1 mM glycine is not a saturating concentration. 
Much better results were obtained when single-channel 
recordings at 10 mM glycine were included in the fits. 
This inclusion of 10 mM of data was essential to obtain 
good agreement of the macroscopic deactivation time 
constant predicted by the fit with our measurements in 
the concentration jump experiments and to predict the 
distinctive, unusually long single-channel activations, 
separated by long shut times. 

For fitting, the idealized records are divided into 
groups of openings that are separated by a shut time 
that is longer than a specified l crll . The aim of this is 
to obtain groups of openings that come from the 



TABLE III 





Empirical fit of mixtures of 


exponential probability density functions to the apparent shut-time distributions 




[Glycine] 


t1 


t2 


t3 


t4 




Area (%) 


Area (%) 


Area (%) 


Area (%) 




ms 


MS 


ms 


ms 


20 yiM n = 2 


0.010 + 0.002 


0.092 + 0.028 


0.589 + 0.027 


6.06 ± 0.01 




(89 ± 3) 


(8±1) 


(3±1) 


(1 ±0.4) 


30 pM n = 2 


0.006 ± 0.001 


0.044 ± 0.014 


0.479 ± 0.055 


5.2 ± 2.5 




(89 ± 2) 


(8 ±2) 


(2 ±0.2) 


(1 ±0.2) 


50 pM n = 2 


0.005 ± 0.003 


0.085 ± 0.071 




2.83 




(99 ± 1) 


(0.9 ±0.8) 




(0.3) 


100 pM n = 3 


0.005 ± 0.001 


0.027 ± 0.007 


0.210 ±0.042 


2.02 ± 0.66 




(98 + 1) 


(2±1) 


(1 ±0.3) 


(0.5 ± 0.2) 


1 mM n = 3 


0.007 + 0.001 


0.045 + 0.014 


0.354 + 0.092 


1.97 ±0.35 




(92 + 2) 


(6 + 1) 


(2 ±0.7) 


(0.9 ± 0.3) 


10 mM n = 3 


0.006 + 0.001 


0.024 + 0.005 


0.219 


1.09 ±0.05 




(85 ± 12) 


(11 ±9) 


(6) 


(1.4±1) 



n is the number of idealized patches used for each concentration. The time constant value and fractional area (in brackets) of each component are the 
average + SDM for each concentration. 
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same individual channel (refer to Materials and 
methods; Colquhoun et al., 1996, 2003). Because 
there was some ambiguity in the choice of t ail (see 
about P op en curve above) , we did several preliminary 
fits with different 4nt values. For most mechanisms 
tested, the fits were consistently more stable when the 
£ cr it was set to 100 ms at 20-30 uM, 10 ms for bursts at 
higher concentrations (up to 1 mM), and 7-10 ms at 
10 mM glycine. 

To test the quality of a model, the predictions from 
a fit were compared with the experimental apparent 
shut-time and apparent open-period distributions. In this 
case, open and shut times were almost independent, 



so conditional distributions of open times were not 
a useful criterion for distinguishing mechanisms 
(Colquhoun and Hawkes, 1987; Colquhoun et al., 
2003; Burzomato et al., 2004). Additional elements 
for evaluating the quality of a mechanism were the 
physical plausibility of the estimated rate constant val- 
ues and their consistency across different sets. Finally, 
simulations of steady-state single-channel records and 
calculations of macroscopic currents (concentration 
jumps) , on the basis of the model and rate constants 
that were fitted, were compared with our experimen- 
tal data. In total, 44 different putative mechanisms 
were fitted (see Discussion). Most failed to describe 
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Figure 3. Fits of Scheme 1 to 0(2 single-channel data. (A) Scheme 1 (Mangin et al., 2003) contains two binding sites and a single open 
state linked to the fully liganded shut state. Directly connected to the bound shut states are short-lived desensitized states. Entry into 
these states is independent of glycine concentration. (B) The results of simultaneous fit of a dataset (containing four single-channel 
patches) with Scheme 1. The fit was performed with the assumption of interacting binding sites. The histograms are the experimental 
open-period and shut-time distributions, the continuous curves are the calculated distributions predicted from the model and fitted 
rate constants (resolution, 20 us), and the dashed lines are the predicted fits at perfect resolution (i.e., when no events are missed). 
(C) Simulation of a single-channel record (at 20 pM glycine) using the fitted rate constants of one of the sets fitted with Scheme 1. This fit 
failed to predict the long shut times between groups of openings we observed in the experimental records, but it instead predicted that 
at low concentration we should observe frequent bursts. (D) Similar to C, using a different dataset (same set as in B), fitted with Scheme 1. 
Although this set can predict long gaps and clusters of openings, it also predicts very short bursts with lower -P open that are not observed 
at any of our single-channel records (compare with Fig. 2 A). 
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the data well, and only a few of the most plausible can- 
didates will be described here. 

Scheme 1: a mechanism with two binding sites and two 
extra shut states 

Kinetic Scheme 1 (Fig. 3 A) has been proposed by 
Mangin et al. (2003) as a minimal kinetic mechanism 
for the a2 homomeric GlyR on the basis of the analysis 
of macroscopic and single-channel currents activated 
by fast agonist application to outside-out patches. 
Scheme 1 contains two short-lived "desensitized" states, 
each connected to a bound shut state. This scheme re- 
sembles closely the one proposed for the activation of 
GABA A receptors byjones and Westbrook (1995) . Good 
fits of apparent open periods and shut times could be 
obtained, but only if the binding sites were allowed 
to interact (i.e., no constraints imposed on the binding 
rates), as for the analogous model for al|3 receptors 
(Burzomato et al., 2004). 

Examples of distributions of apparent open periods 
and shut times for a set fitted with Scheme 1 are plotted in 
Fig. 3 B. As shown from this example, the fit of Scheme 1 
predicts well the shape of the apparent shut- and open- 
time distributions. The results from the fits to three data- 
sets are summarized in Table IV. As seen from the table, 



TABLE IV 

Estimated parameters from Scheme 1 fitted to a 2 GlyR single-channel 
data 



Parameter name and unit 


Scheme 1 


Rate constants ± CVM (%) 




a (s- 1 ) 


330 ± 17% 




120,000 ± 4% 




3,450 + 53% 


A- (s- 1 ) 


176 ±71% 




3,860 ± 10% 




2,740 ± 20% 


hi (M _1 s _1 ) 


18,200 ± 59% 


k-i (s- 1 ) 


1,480 ± 49% 


h-2 (M _1 s _1 ) 


4.1 x 10 s ± 75% 


k- s Or 1 ) 


6,680 ± 49% 


Equilibrium constants ± CVM (%) 




E = p/a 


384 ± 13% 


ft = d u / d\- 


28 ± 63% 


ft = d 2t / ck- 


1.64 ±36% 


K, = k-i/hi (M-M) 


107 ± 73% 


K 2 = k- S /hi (uM) 


571 ± 97% 


P open curve parameters (±SDM) predicted from this scheme 


max P open 


0.995 + 0.001 


EQo (uM) 


134 ±55 


»ii 


1.94 ±0.03 



Rate constants are means of fits to three data sets ± CVM. The equilibrium 
constants (E, ft, ft, Ki, and K?) are the means of the three values calculated 
for each separate set. Also listed are the average max P open , £C 50 , and Hill 
slope calculated from the fitted rate constants (±SDM). These fits have 10 
free parameters. 



the values for the gating rate constants (the shutting rate 
constant a and the opening rate constant (3) were highly 
consistent across the three sets, with a very low coeffi- 
cient of variation. 

Rate constants other than those for opening and shut- 
ting were highly variable across the three sets (coeffi- 
cient of variation [CVM] of 50% or more), suggesting 
that the scheme cannot estimate these rate constants 
reliably. Note also that the ECsq predicted from the fit- 
ted rate constants was almost sixfold lower than the EC^ 
measured in outside-out recordings. To test further the 
correctness of Scheme 1 , we calculated the macroscopic 
current expected in response to 2-ms fast concentration 
jumps with 10 mM glycine. The current response was 
calculated for each set of fitted rate constants. The pre- 
dicted rise time of macroscopic currents from the three 
sets fitted had a mean time constant of 10 ± 5 ms, which is 
^25-fold slower than the 20-80% rise time measured from 
currents evoked by 10 mM glycine (t = 0.42 ± 0.08 ms; 
Table I) . On the other hand, the mean deactivation 
time constant predicted by Scheme 1 was 120 ± 70 ms 
(n = 3 sets), which is very close to the experimental value 
of 130 ± 10 ms (the slow component; Table I). This 
corresponds to a mean burst length calculated from 
Scheme 1 of 120 ± 70 ms at 1 nM glycine (n = 3 sets). 

We also checked whether simulated single-channel 
records (SCSIM program) resembled the pattern seen 
in experimental records. Single-channel simulations at 
20 uM glycine using Scheme 1 failed to reproduce the 
most striking features we observed in our single-channel 
experimental records, namely the presence of groups 
of openings with duration of up to a few seconds long, 
separated by very long gaps lasting for minutes. For two 
sets, bursts (with lengths ranging from tens of millisec- 
onds to ^1 s) appeared continuously in the simulated 
records, and there was no sign of the very long shut 
times between them, as seen in experimental records 
(compare Fig. 3 C with Fig. 2 A) . Only the simulation from 
the third set showed long gaps (Fig. 3 D), but even this 
realization did not show long, high P ope n groups of open- 
ings. The mean burst length calculated from Scheme 1 
at 20 pM glycine was only 300 ± 30 ms (n = 3 sets). These 
simulations show that Scheme 1 could not describe our 
experimental results adequately, despite a good 
match to the dwell-time distributions. Therefore, we 
proceeded to test alternative mechanisms for our data. 

Scheme 2: a flip mechanism with two binding sites and one 
open state 

Kinetic Scheme 2 (Fig. 4 A) is a variant of the flip mech- 
anism, originally proposed by Burzomato et al. (2004) 
for al(3 and al GlyRs. Scheme 2 is a modified version 
containing two binding sites (rather than the three needed 
to reach maximum P opl , n in al homomers) . We have also 
chosen to omit from Scheme 2 the two partially liganded 
open states included in the al/al(3 flip mechanism. 
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Including more than one open state in the mechanism 
barely improved the quality of the fits to the present 
data and resulted in poorly determined rate constant 
estimates for the entry and exit from partially liganded 
open states. This probably means that openings from 
partially liganded receptors are rare, so there are not 
enough transitions in our records to be able to describe 
these states accurately. The predominance of visits to a 
single open state is confirmed by the lack of any substan- 
tial correlation between open and shut times, and by the 
fact that the distribution of apparent open-period dura- 
tions is not far from a single exponential. 

An example of a dataset fit with Scheme 2 is shown in 
Fig. 4 B. There is a good agreement between the experi- 
mental distributions (histograms) and those that are 



predicted from the fits using the model (continuous 
line). The set shown in Fig. 4 B is the same as in Fig. 3 B 
for comparison. The average rate constants from three 
different sets are summarized in Table V. Note that 
these fits were done with the assumption that binding 
sites are independent; i.e., they do not interact, so all 
binding steps to a given conformation (R or F) are iden- 
tical, regardless of how many sites are already occupied. 
This was done by imposing the constraint that k+i = k+? 
and k.\ = k_9. 

As seen from the values in Table V, the fits with 
Scheme 2 are very consistent across the three sets, and 
most rate constants had low coefficients of variation. 
The least well-defined rate constants are those for enter- 
ing and exiting the singly occupied flip state (Si and y{) , 
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Figure 4. Fits with a flip model containing two binding sites (Scheme 2) . (A) Scheme 2 is a flip model with two binding sites and a single 
open state that can be accessed from the fully liganded shut state after flipping. (B) Plots that show open-period and shut-time distribu- 
tions fitted with Scheme 2. The same set of four patches is used as in Fig. 3 B for comparison. The fits with Scheme 2 are performed 
using constraints for non-interacting binding sites (resolution, 20 ps) . (C and D ) Simulation of a single-channel record (at 20 pM glycine) 
using the fitted rate constants of two of the sets fitted with Scheme 2 (the set in D is the one shown in the histograms). The predictions 
from all sets fitted with this model are in a very good agreement with the observed single-channel records (compare with Fig. 2 A). 
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although they are clearly quite slow. This is not surpris- 
ing given the evidence already cited that partially ligan- 
ded states occur only rarely. Scheme 2 again predicts 
that the efficacy for the opening step is very high (350) , 
and that the efficacy for the flipping steps increases 
as the number of bound glycine molecules increases 
(Table V) . Also, the shut flipped states have a 40-60-fold 
higher affinity for glycine than the shut resting states. 
This is a well-characterized feature of the flip mecha- 
nism, which suggests that the binding of the agonist sta- 
bilizes the receptor in a higher affinity form reached by 
a conformational change that precedes channel open- 
ing (Burzomato et al., 2004). The average maximum 
^open predicted by the fits to Scheme 2 was 0.99, and the 
predicted nn was 1.89, in good agreement with the value 
of 1.99 (ra H ; Fig. 1 D) obtained from our macroscopic 
dose-response curve experiments (see also Mangin 
et al., 2003) and with the maximum P ope n of 0.99 from 
single-channel records. Note, however, that there was 
about a threefold difference between the experimental 
macroscopic -EC50 (770 uM) and the one predicted by 
the model (240 uM; Table V) . 



Using the rate constants obtained from the fits de- 
scribed above, we calculated macroscopic current re- 
sponses to a 2-ms fast concentration jump with 10 mM 
glycine. Scheme 2 predicted a time constant of 171 ±6 ms 
(n = 3 sets) for the decay of these currents, which is 
similar to the 130 ± 10 ms we obtained from our mac- 
roscopic currents in outside-out patches. Scheme 2 pre- 
dicted two significant components for the activation 
phase of the current macroscopic current (1.6 ± 0.2 
and 0.19 ± 0.02 ms), but the contribution of the faster 
component was minimal. Consequently, the weighted 
calculated rise time (t = 1.5 + 0.2 ms) reflected mosdy 
the slow component and was almost fourfold slower 
than the rise time measured from currents evoked 
by 10 mM glycine (20-80% rise time, 0.42 + 0.08 ms; 
Table I). 

Finally, steady-state single-channel records (at 20 uM 
glycine) simulated using the fitted rate constants showed 
groups of openings with a mean burst length of 490 ± 
70 ms, separated by long shut times (>1 min; Fig. 4, 
C and D) . These features of single channels are similar 
to our cell-attached recordings. It is not surprising that we 



TABLE V 

Estimated parameters from Schemes 2, 3a, and 3b fitted to a 2 GlyR single-channel data 



Parameter name and unit 



Scheme 2 



Scheme 3a 



Scheme 3b 



Rate constants ± CVM (%) 

a (s- 1 ) 

7 3 (s" 1 ) 
8s(0 
V (M _1 s _1 ) 
h- (s" 1 ) 
K (M -1 s -1 ) 
k- (s- 1 ) 

Equilibrium constants ± CVM (%) 

E=j3/a 
Fi = 8 l/ 7 1 
F 2 = S 2 / y 2 
Fs = 8 3/ y 3 
K R =k-/h (uM) 
Kg = (uM) 
Kr/K ¥ = F,/F 1 

Predicted P opcn curve parameters ± SDM 

max P open 
EQo (uM) 
»ii 



366+ 12% 
127,000 + 9% 
4 ±81% 
0.134 ± 76% 
4,418 + 11% 
4,316 ± 12% 



7.9 x 10 6 ± 27% 
677 ± 12% 
1.04 x 10 5 ±5% 
439 ± 22% 

347 ± 6% 
0.031 ± 28% 
0.977 ± 24% 

4,241 ± 26% 
86 ± 37% 
46.5 ± 28% 

0.990 ± 0.006 
243 ± 52 
1.89 ± 0.005 



320 + 3% 
119,919 ±7% 
292 ± 98% 
0.014 + 82% 

2.1 x 10"' ±48% 
1.7 x 10~ 9 ±42% 

4,109 ± 20% 
2,572 ± 38% 
3.1 x 10 7 ± 61% 
719 ± 37% 
0.779 x 10 5 ± 21% 
154 ± 36% 

373 ± 4% 

2.2 x 10"" ±63% 
0.011 ±47% 
0.573 ±31% 
1,879 ±21% 

31 ±21% 



0.993 ± 0.002 

389 ± 66 
2.49 ± 0.06 



321 + 4% 
120,000 ± 8% 



106 ± 100% 
0.36 ± 100% 
4,152 ± 19% 
2,605 + 38% 
2.4 x 10' ± 75% 
659 ± 42% 
0.781 x 10 5 ± 26% 
158 ± 38% 

372 ± 4% 

0.017 ± 52% 
0.570 ± 29% 
1,919 ± 20% 
48 ± 29% 



0.993 ± 0.002 
400 ± 67 
2.49 ± 0.05 



Rate constants are means of fits to three data sets ± CVM. The equilibrium constants (E, Fi, F%, F$, K R , and K^) are calculated for each separate fit, and 
means are listed here (because of the variability, the means of the equilibrium constants are not exactly the same as the ratios of the mean rate constants). 
Likewise, the mean values of the relative affinity of glycine for resting and flipped (K R / K F ) are given. Also listed are the average max P open , -EC50, and 
Hill slope calculated from the fitted rate constants (±SDM). The number of free parameters: 9 out of 12 in Scheme 2, 10 out of 18 in Scheme 3a, and 9 
out of 14 in Scheme 3b. Other rate constants are set by microscopic reversibility or constrained by the assumption that the binding sites are equivalent 
and independent. 
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saw no channel activity at 10 uM glycine, because cal- 
culations based on the fitted rate constants predict 
that activations should be separated by long silent peri- 
ods of >4 min (see Discussion). 

To summarize, based on the above observations and 
on preliminary fits to different models, Scheme 2 is the 
minimal model that can describe best both our single- 
channel and macroscopic data. 

Scheme 3: a flip mechanism with three binding sites 
and one open state 

The mechanism in Scheme 2 contains two binding sites. 
This is different from what we suggested earlier for the 
ctl homomeric GlyR, namely that maximal activation of 
the receptor requires the occupancy of three binding 



sites by glycine (Beato et al., 2004). It is of course pos- 
sible that the number of functional binding sites is 
different between a2 and al homomers, but this is an 
important point to test. We therefore extended Scheme 2 
and tested two related mechanisms (Scheme 3a and 
Scheme 3b), shown in Fig. 5 A. The two mechanisms 
differ only in that the monoliganded flipped state AF is 
omitted from Scheme 3b. As for Scheme 2, we assumed 
that binding sites are identical (within each form of the 
receptor) and do not interact. This was imposed by 
using the following constraints: k+i = k + 9 = k+ 3 and k_i = 
k-z = &_3 for both R and F states. Note that because of 
these constraints, the number of free parameters in 
Schemes 2 and 3b is the same, and both have nine free 
rate constants. 
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Figure 5. Fits with flip models containing three binding sites (Schemes 3a and 3b) . (A) Scheme 3a (top) differs from Scheme 2 (Fig. 4 A) 
in the presence of a third agonist binding site. Scheme 3b (bottom) is a variant that lacks the partially liganded AF state. Both models 
have only one open state that can be accessed after flipping, when the three binding sites are occupied by glycine. (B) Plots that show 
open-period and shut-time distributions fitted with Scheme 3a. Identical fits were obtained with Scheme 3b. Note also the similarity with 
predictions from Scheme 2 (Fig. 4 B) . Fits with models 3a and 3b were performed with constraints for non-interacting binding sites 
(resolution, 20 us). (C) Simulation of a single-channel record (at 20 pM glycine) using the fitted rate constants of the set shown in B. 
Similar simulations are obtained for Scheme 3b. The prediction from both models is similar to those from Scheme 2 and in a very good 
agreement with our cell-attached data (compare with Figs. 2 A and 4, C and D) . 
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The apparent quality of the predictions of the dwell- 
time distributions was good and very similar for the two 
models; therefore, we show only the fits of Scheme 3a in 
Fig. 5 B, using the same set of data as in the previous fig- 
ures for comparison. The predictions of the dwell-time 
distributions were almost identical for Schemes 2, 3a, 
and 3b (compare Fig. 4 B and 5 B), and the values of 
the rate constants obtained by fitting the three schemes 
(see Table V) were very similar, although the coefficient 
of variation was overall higher for the fits with Scheme 3b. 
The most variable rate constants are again those that de- 
scribe the transitions between partially liganded flipped 
and resting states. 

Scheme 2 is favored over Schemes 3a and 3b because 
of its better prediction of the Hill slope (1.89 ± 0.01 for 
Scheme 2, but 2.49 ± 0.06 for Scheme 3a and 2.49 ± 0.05 
for Scheme 3b; see also Table V) . Because we were not 
able to obtain a one-channel P op en curve, our only ex- 
perimental estimate of a slope (wh = 1.99; Fig. 1 D) 
comes from the concentration-response curve obtained 
from the macroscopic current records in the outside- 
out configuration. We do not know how much this mea- 
surement is distorted by desensitization, so this is not a 
conclusive argument. 

To check further the quality of Schemes 3a and 3b, we 
calculated macroscopic currents and simulated single- 
channel currents. Fig. 5 C shows a simulated single- 
channel trace (20 uM glycine), using the rate constants 
from a set fitted with Scheme 3a. In this record, open- 
ings occur in groups several hundred milliseconds long, 
separated by very long shut times, as in our experimen- 
tal traces. Similar results were obtained from simula- 
tions with Scheme 3b (not depicted). Therefore, as with 
Scheme 2, Schemes 3a and 3b predicted well our single- 
channel data. The rise and deactivation time constants 
of macroscopic currents were calculated from the fitted 
rate constants. The average weighted rise time of a cal- 
culated current activated by 2 ms of 10 mM glycine was 
1.5 ± 0.4 ms (n = 3 sets) , and the deactivation time con- 
stant was 150 ± 30 ms. Both time constants were almost 
indistinguishable from those calculated for Scheme 2. 
The relationship between burst length and macroscopic 
deactivation is considered in the Discussion. 

DISCUSSION 

The main aim of this paper is to characterize the activa- 
tion mechanism of homomeric a2 GlyRs. We used the 
method of maximum likelihood fitting of kinetic mech- 
anisms to idealized single-channel data (Colquhoun 
et al., 2003) for the a2 GlyR. In brief: (a) single-channel 
data are idealized by time-course fitting, and patches of 
different agonist concentrations are grouped into sets; 
(b) a kinetic mechanism is postulated (usually taking 
into account exponential fits to experimental dwell-time 
distributions) and initial guesses are given for the rate 



constants; (c) the idealized records, the appropriate re- 
solution, and the postulated mechanism are introduced 
into a program (HJCFIT) that can calculate the likeli- 
hood (i.e., the probability of the entire sequence of ap- 
parent open and shut times, given the current estimates 
of the rate constants) and maximize it by optimizing 
rate constants; and finally, (d) the validity of the model 
and of the final rate constants is checked by comparing 
their predictions with the experimental data (dwell-time 
histograms). For the a2 GlyR, we also examined how 
well the mechanism and rates fitted to single channels 
predicted the time course and concentration dependence 
of the macroscopic currents elicited by concentration 
jumps. The terms activation, burst, and cluster are de- 
fined at the beginning of the Results section. 

Properties of a2 GlyR single-channel activity 

The first kinetic analysis on GlyRs, performed on cul- 
tures of spinal motorneurons, showed a characteristic 
feature of these channels: the frequency of long and 
short openings is concentration dependent. In record- 
ings at low glycine, increasing glycine concentration 
increases the proportion of long openings (without 
changing the time constants of the different opening 
components; Twyman and Macdonald, 1991). The same 
phenomenon was subsequently detected in recombi- 
nant homomeric and heteromeric GlyRs that contain 
the al subunit, where bursts at higher glycine concen- 
trations are longer because they contain more openings 
and because the openings themselves appear to be 
longer (Fucile et al., 1999; Beato et al., 2002, 2004; 
Burzomato et al., 2004). 

This is not the case for the activity of a2 GlyRs, where 
openings occur in long groups, regardless of concentra- 
tion (Fig. 2 A) . At first, the length of the groups of open- 
ings (Fig. 2 A) misled us into thinking that they were 
clusters (as defined above). Even at the lowest concen- 
tration of glycine at which channel activity could be ob- 
tained (20 uM), long groups of openings were observed 
(several hundreds of milliseconds) , and the probability 
of being open within the group of openings was high 
(0.96; tent, 100 ms). However, the macroscopic concen- 
tration-response curve (Fig. 1 D) suggests that 20 uM 
glycine should produce only a very small response and 
therefore a very low overall single-channel P op en (<5%) . 
This observation, together with the fact that we could 
not observe activity at concentrations lower than 20 pM, 
implies that the long groups of openings at 20 pM 
glycine must be (unusually long) individual activations 
of the channel (possibly lengthened a bit by rebinding) . 
In other words, the true P op en must be low at 20 pM 
and should be measured including the very long shut 
periods that separate the groups of openings. This 
cannot be done because there is no way to tell whether 
two adjacent activations originate from the same indi- 
vidual channel. 
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The use of concentration jumps for fast agonist ap- 
plication provides the benefit of mimicking what oc- 
curs in vivo: a fast and brief exposure of the receptors 
to saturating concentrations of neurotransmitter in 
the cleft during synaptic transmission (Clements et al., 
1992; Colquhoun et al., 1992) . Moreover, when the ap- 
plication of agonist is short (1-2 ms), the decay of the 
current is known to reflect the duration of the burst at 
very low agonist concentration, i.e., the duration of 
channel activations (Wyllie et al., 1998). When jumps 
of 2 ms to three different glycine concentrations were 
applied, two deactivation components were identified. 
The time constant of the slow component, which is by 
far the predominant one, was ^130 ms, which is ^17 
times slower than the deactivation of al(3 channels 
(7.5 ms with 10 mM intracellular CP; Pitt et al., 2008) . 
Note that Mangin et al. (2003) reported for a2 cur- 
rents, using high chloride in the pipette, a deactiva- 
tion time course similar to the one we observed (153 ms) . 
We found (unpublished data) that changing intra- 
cellular chloride concentration to 130 mM had litde 
effect on a2 currents, in contrast with our findings 
for oil|3. 

Choice of mechanism 

Despite the fact that the single-channel records look 
very different from those seen with al-containing ho- 
momers or heteromers, it was found eventually that the 
best description among the 44 mechanisms that we 
tested was provided by the flip model (Burzomato et al., 
2004; Lape et al., 2008). An overview of the type of 
mechanisms tested is shown in Fig. 6. 

The fit of the shut-time distributions was not perfect 
with this mechanism, and neither was the prediction of 
the rise time for concentration jumps. However, both 
were better than could be obtained with the main rival 
class of mechanisms, as represented here by the Mangin/ 
Jones-Westbrook type of model (Scheme 1). In our 
original work on the flip model (Burzomato et al., 
2004) , the Jones-Westbrook-type mechanism fitted at 
least as well as flip, albeit with four more free parame- 
ters. The preference for the flip model then was based 
to a certain extent on physical grounds (it allowed a fit 
without having to postulate interaction between differ- 
ent binding sites) , partly because the flip model had a 
plausible structural interpretation (C-loop capping, since 
confirmed in muscle nicotinic channels by Mukhtasimova 
et al., 2009). Mechanisms with preopening intermedi- 
ate states work well for wild-type and mutant GlyR and 
muscle nicotinic channels and for GABA A channels 
(Burzomato et al., 2004; Plested et al., 2007; Lape et al., 
2008; Mukhtasimova et al., 2009; Keramidas and Harrison, 
2010). In the present case, there was more reason to 
discard a Mangin/Jones-Westbrook type of model in 
favor of flip than in our earlier work. For example, the 
fits with flip give more consistent rate constant values 



across experimental sets. Flip is also somewhat better at 
predicting EC50 values, the rise time of macroscopic 
currents, and the pattern of channel activity in single- 
channel records. 

There are two respects in which the fit with the flip 
mechanism (Scheme 2), despite being the best we 
found, is less than perfect. The EC50 predicted by the fit 
to Scheme 2 was almost three times smaller than that 
observed for the macroscopic concentration-response 
curve. Also, the predicted rise time for responses to 
pulses of 10 mM glycine was over three times slower 
than the observed value (although that is a great deal 
closer than for Scheme 1, for which it was 25 times too 
slow) . It could well be that the flip model is too simple 
to account for these details. The least plausible aspect 
of the flip mechanism is the postulate that flipping is 
a concerted conformation change. It might be better 
to allow flipping (or "priming") to occur separately in 
each subunit (Mukhtasimova et al., 2009). Desirable 
though that might be, in its general form, the primed 
model has so many free parameters that they cannot 
all be estimated from the data. We have preferred an 
approximation that allows all of the rate constants to 
be estimated. 

Features of the activation mechanism for ct2 GlyRs: 
the number of open states and binding sites 
The general flip model that describes the activation of 
homomeric al and heteromeric alp GlyRs contains 
three open states (Burzomato et al., 2004). This is in 
line with early data on native and recombinant recep- 
tors that detected at least three distinguishable open 
states (Twyman and Macdonald, 1991; Beato et al., 
2002, 2004; Lewis et al., 2003). In contrast to the al- 
containing channels, a2 receptors could be fitted quite 
well with a flip model that contains only one open 
state. This suggests that partially liganded openings are 
either rare or too short to be detected. The latter expla- 
nation is supported by preliminary fits with a version of 
Scheme 2 that included monoliganded openings (not 
depicted) . Partially liganded openings had a very high 
closing rate constant and were so short that, if they truly 
occurred, most of them would be missed with our reso- 
lution of 20 or 25 ps. For these reasons, we considered 
it more sensible to leave only one fully liganded open 
state. Another reason for this choice is that we could not 
detect open-shut correlations in our experiments with 
a2 receptors, and these are not expected to occur in 
mechanisms with only one open state (Colquhoun and 
Hawkes, 1987). 

In the general flip scheme used to fit homomeric 
al channels, maximum P ope n is reached when three 
(out of the possible five) binding sites are occupied by 
glycine. For the a2 channel, good, almost identical 
fits of single-channel data were obtained with flip 
mechanisms with two or three binding sites, and these 
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resulted in very similar values for the rate constants. 
The biggest difference between the predictions of 
mechanisms with two or three binding sites was in the 
value of the Hill slope of the concentration-response 
curve. Scheme 2 (two binding sites) predicted a value 
of 1.9, in better agreement with our experimental 
data from macroscopic currents in the outside-out 
configuration (1.99; Fig. 1 D). The limitation is that 
we only have a macroscopic concentration-response 
curve, which is less satisfactory than a single-channel 
P op en curve because it is susceptible to errors from de- 
sensitization, as well as possible differences between 
outside-out patches and cell-attached patches. Even 
so, it is clear that occupancy of all the five potential 
agonist binding sites in a homomeric nicotinic-type 
receptor is not needed to achieve the maximum P op en, 
in agreement with our own conclusions for the ctl 
GlyR homomer (Beato et al., 2004) and with those of the 
elegant experiments of Rayes et al. (2009) on a7-5HT 3 
chimeric receptors. 



Features of the activation mechanism of a2 GlyRs: 
efficacy and affinity 

For the purpose of characterizing the kinetic behavior 
of homomeric a2 GlyRs, steady-state single-channel data 
were fitted to 44 different putative reaction schemes, 
only four of which are described here. These potential 
mechanisms fall into two distinct types. In Scheme 1, 
originally proposed by Mangin et al. (2003), the satu- 
rated receptor has only two states, open and shut, and 
the efficacy of the agonist depends only on the equilib- 
rium constant for this step, E = fti/oL-z, in the classical 
way. The maximum P ope n is 



P open (max) = — — . 



(1) 



To get good fits, it was necessary to add two extra shut 
states, in the manner of Jones and Westbrook (1995). 
Although these states are labeled as "desensitized," they 
are quite short-lived, and there is no reason to think 
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Figure 6. Activation mechanisms tested by global fitting to single-channel data. The mechanisms shown here are the most general 
forms of each of the 44 mechanisms we tested. In all schemes, the letter A denotes an agonist molecule and its subscript indicates the 
number of agonist molecules bound to the receptor. The letters R and R* denote resting shut states and open states of the receptor, 
respectively. (A) Sequential scheme with five binding steps and open states (nine "subset" schemes tested in all, with two to five binding 
steps, allowing opening either from each bound state or only from the one or two highest liganded states) . (B) Jones and Westbrook- 
type scheme (Jones and Westbrook, 1995), with an open state and a distal desensitized state (indicated by D) for each binding step (five 
schemes tested with one to three binding steps, each leading to a desensitized state, allowing the channel to open only from the highest 
liganded states). (C) Flip-type scheme, where an intermediate flipped shut state (F) connects resting and open states (Burzomato et al., 
2004; five schemes tested with one to five binding steps, a single flipped state, and a single open state; eight schemes tested, where the 
number of flipped and open states was either the same as the number of binding steps or smaller; three schemes tested, where, after the 
first binding step, the agonist can bind only to flipped states; three schemes with additional flipped states). (D) A subset of the primed 
mechanism (generalization of the flipped mechanism [Mukhtasimova et al., 2009]; 11 schemes tested with two to three binding steps 
and fewer open states than the maximum shown here) . Several independent priming steps can occur. Here, the channel can open after 
one, two, or three priming steps (indicated as F, FF, or FFF) . 
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that they have anything to do with macroscopic desensi- 
tization. The extra shut states are, therefore, entirely 
arbitrary. Furthermore, to obtain a good fit, it was nec- 
essary to assume that the two binding sites in Scheme 1 
interact (the binding constant for the second binding is 
not the same as for the first) . 

For both of these reasons, we turned to flip mecha- 
nisms (Schemes 2, 3a, and 3b) . In this class of mecha- 
nisms, the binding sites are supposed to be independent, 
and there is an extra shut state (the flipped conforma- 
tion) between resting and open states. Therefore, the 
fully saturated receptor has three states (resting, flipped, 
and open) rather than two, and the maximum P ope n de- 
pends on two different equilibrium constants. These 
are the equilibrium constant, E, for the open-shut reac- 
tion of the fully liganded receptor (the "opening effi- 
cacy"), and the equilibrium constant for the preceding 
conformation change from resting to flipped state, F 
(the "flipping efficacy"), 



^openCmax) 



EE 



l + F{\ + E) 



(2) 



This can be written in the classical form if we define an 
"effective efficacy," E eB , as 



E eS =E 



l + F 



(3) 



With this definition, the maximum P open for the flip type 
mechanism can be written in the form of Eq. 1 as: 



^open (max) 



(4) 



The flip mechanism was first proposed for glycine 
al(3 receptors by Burzomato et al. (2004). It follows 
from Eq. 3 that if the flipping efficacy, F, is small, the 
overall efficacy may be much lower than expected from 
the opening efficacy, E, and Lape et al. (2008) gave evi- 
dence that a reluctance to flip (small value of F) was the 
basis of partial agonism for both glycine and nicotinic 
receptors. A similar explanation for the partial agonism 
of choline at nicotinic receptors was proposed by Lape 
etal. (2009). 

The fits of the flip mechanism (Scheme 2) suggested 
that glycine has a 46-fold higher affinity for the flipped 
conformation than for the resting conformation, little 
different from the 65-fold higher affinity for al(3 GlyRs 
found by Burzomato et al. (2004). The results resemble 
al(3 GlyRs also in that the increase in affinity is brought 
about almost entirely by an increase in the association 
rate constant for the flipped conformation, rather than 
a decrease in the dissociation rate. The association rate 
constant for the resting conformation, ^10 5 M _1 s _1 , is 
unusually low. Even the association rate constant for the 
flipped conformation, ^8 x 10 6 M _1 s _1 , is still slower 



than that for the binding of acetylcholine to muscle- 
type nicotinic receptors, which is >10 8 M _1 s _1 . Although 
the high affinity ratio (between resting and flipped 
states) results in a high degree of apparent cooperativ- 
ity (the Hill slope is close to its maximum value of 2), 
the absolute value of the flipping equilibrium constant 
is only near 1, compared with 27 for a 1(3 GlyRs (Lape 
et al., 2008). Nevertheless, the overall efficacy is high 
because of the high equilibrium constant (E>) for the 
shut-open step. 

The values of the gating rate constants were the most 
consistent across experimental sets in our maximum 
likelihood fitting, regardless of which model was used 
for fitting. They ranged between 300 and 400 s" 1 for the 
shutting rate a and 110,000 and 130,000 s 1 for the 
opening rate (3. The values given in Table V for the fit of 
Scheme 2 have the following implications. The fast 
opening rate constant is very similar to that for the 
heteromeric al(3 receptors (Burzomato et al., 2004), 
but the shutting rate constant is 20 times slower, so the 
mean length of the individual openings is ^3 rather 
than 0.14 ms for heteromeric al(3 receptors. The mean 
efficacy for the opening step is ^350, much higher 
than the gating efficacy of either al(3 heteromeric GlyR 
(a = 7,000 s" 1 ; (3 = 129,000 s" 1 ; E= 20) or al homomeric 
GlyR (a = 700 s -1 ; (3 = 28,000 s _1 ; E = 38; Burzomato 
et al., 2004). In al-containing receptors, the lower (but 
still large) efficacy for the open-shut step stems from 
a faster shutting rate and implies a short mean open 
time for the al heteromer, but in the al homomer, it 
appears to result from a slower opening rate constant. 
For a2 receptors, the flipping equilibrium constant for 
saturated receptors (Table V) is around 1, so shut satu- 
rated receptors spend about half their time in the rest- 
ing conformation (A 2 R) and half the time flipped (A 2 F) . 
The overall effective efficacy is still very high, E eff =171, 
from Eq. 3, which implies a maximum P ope n of 99.4%, 
from Eq. 4. 

The energy diagram in Fig. 7 plots the energy barriers 
for the transitions of fully liganded receptors from 
resting to flipped to open for the GlyRs we have so far 
analyzed (al, al (3, and a2) . It is interesting to note that 
the height of the energy barrier that the channel has to 
cross to go from resting to flipped seems to be linked in 
our series to whether the principal subunit in the recep- 
tor is al or a2. It is plausible that the subunit that pro- 
vides the principal side of the binding site should play a 
major part in this transition to flipped, a transition that 
is likely to correspond to the conformational changes 
that follow agonist binding and spread through the ex- 
tracellular domain to reach the channel gate (Grosman 
et al., 2000; Lee and Sine, 2005). Almost all of the se- 
quence differences between al and a2 are in their ex- 
tracellular domain, whereas there is only one amino 
acid difference in M2, the main domain that lines the 
channel and gates it shut. On the other hand, the M2 
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sequence is less conserved in the (3 subunit. Notably, 
the energy difference between the initial resting state 
and the final open state is very similar for the three 
receptors in the series, corresponding to the fact that 
the effective efficacy E eS (Eq. 3) of glycine is very high 
for all three. 

The mechanistic interpretation of a2 GlyR bursts 

It was shown rigorously by Wyllie et al. (1998; appendix) 
that the time constants (and mean length) for bursts 
(activations) measured at sufficiendy low agonist con- 
centrations must be the same as the time constants for 
the decay of macroscopic concentration jumps to zero 
concentration. Obviously, bursts cannot be measured at 
zero concentration, so this theorem begs the question, 
What concentration is "sufficiently low"? For most sorts 
of receptor, it is possible to record bursts at a concentra- 
tion that is low enough for there to be good agreement 
between the single-channel mean burst length and the 
macroscopic deactivation time constant. According to 
our interpretation, the a2 receptor is an exception to 
this rule. 

The observed time constant for the predominant slow 
component of the macroscopic deactivation was 130 ms 
after a step from 10 mM glycine to zero (Table I). The 
fit of Scheme 2 predicts a predominant deactivation 
time constant of 171 ms, not far from that observed. On 
the other hand, this value does not agree very well with 
the observed mean burst length of 490 ms measured at 
the lowest concentration of glycine at which measure- 
ments could be made, 20 pM. Evidendy, quite a lot of 
reopening from monoliganded channels must occur at 
20 pM glycine, and this reopening prolongs the burst 
length to a value that is well above its low concentration 
limit, even at 20 pM glycine. Calculations from the fitted 
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Figure 7. Energy diagram for the resting-flipped-open transi- 
tions for the fully liganded al, oil(3, and a2 GlyR. Calculations 
used a frequency factor of 10~ 7 s _1 (Andersen, 1999) and the re- 
sults of the fit of Scheme 2 (flip with two binding sites and one 
open state) for a2 GlyR, and data from Burzomato et al. (2004) 
for the al-containing GlyRs. 



rate constants for Scheme 2 predicts that at a really low 
glycine concentration, 1 nM, the mean burst length 
should be 171 ms, in agreement with the predicted de- 
activation rate for Scheme 2. Clearly, the way the num- 
bers work out in this case makes it impossible to record 
bursts at a sufficiently low concentration. 

To make sense of the numbers, it is useful to define 
the relative probabilities that a channel moves next to 
each of the adjacent states in the scheme, regardless of 
how long it takes to do so. For example, in Scheme 2 
(Fig. 8), a monoliganded resting channel AR (state 5) 
can move next to the resting state, R (state 6), it can flip 
to AF (state 3) , or it can bind another glycine molecule 
to get to A 2 R (state 4) . The probabilities are simply pro- 
portional to the rate constants for leaving state 5 
(Colquhoun and Hawkes, 1995; p. 415), so, for exam- 
ple, the probability that a monoliganded channel disso- 
ciates is 



(5) 



where c is the concentration of glycine. 

The concentration dependence of the length of an 
activation arises, in the case of ot2 receptors, largely 
from the binding of glycine to the monoliganded flip 
state (AF). When 20 pM glycine is applied, opening will 
take place usually via A 2 R because the monoliganded 
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Figure 8. Features of Scheme 2 as fitted to a2 GlyR data. The 
diagram shows the flip mechanism with two binding sites fit- 
ted, with the different states numbered to provide a key to the 
physical interpretation of the results in the Discussion. Numbers 
above the arrows indicating the transitions are the average rate 
constants from the fits. Units are s _1 or M - ^ -1 , as appropriate 
(see Table V) . 
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flipping rate is small (ir 54 tt 42 is almost 20 times greater 
than tt 53 tt 32 ) . After an opening, the channel must re- 
turn to A 2 F (state 2), from where it is likely to reopen 
(ir 2 i = 0.96), but if it does not reopen, it is three times 
less likely to go to AF (state 3) than to unflip to A 2 R. 

If this occurs, even at the lowest concentration of 
glycine for which recordings could be made, 20 uM, the 
probability that from AF the channel can bind glycine 
again is high (tt 32 = 0.98). Because this binding returns 
the channel to A 2 F, and from there the channel is very 
likely to (re) open, this path accounts for the lengthening 
of the activation at 20 uM glycine. In contrast to that, the 
probability of a second binding occurring to the resting 
conformation, AR, is very low (ir 54 = 0.004 in 20 uM gly- 
cine), and the probability of flipping is even lower. 

There are two reasons for the unusually long activa- 
tions produced by the a.2 receptor even at the lowest 
concentrations. The first is that individual openings are 
longer, on average, for the a2 receptor (3 ms; Table V) 
than for the al(3 receptor (0.14 ms; Lape et al., 2008). 
The second reason is the large number of openings per 
activation seen with the a2 receptor, around 100, even 
at very low concentrations. 

Shut times within an activation turn out to have a 
rather simple interpretation in the case of the a2 recep- 
tor. The distribution of shut times calculated from the 
fitted rate constants for Scheme 2 (Table V) , with Eq. 3.90 
from Colquhoun and Hawkes (1982), is dominated at 
all concentrations by a very short component with a 
mean of ^8 us, close to the mean duration of a sojourn 
in A 2 F (state 2) . Clearly, most of the time during an acti- 
vation is spent oscillating between A 2 F* and A 2 F. There 
is also a small component (a few percent in area) of 
shut times within activations that has a time constant 
that is close to the mean lifetime in A 2 R (state 4) , M).2 ms. 
This clearly represents channels that undergo transi- 
tions A 2 F*-»A 2 F-»A 2 R->A 2 F->A 2 F*. There will very 
rarely be oscillations between A 2 F and A 2 R, so this time 
constant of the shut-time distribution is close to the 
time spent in A 2 R (plus two brief sojourns in A 2 F) . Still 
more rarely (at 20 uM glycine), there will be transitions 
A 2 F* ->A 2 F->AF->A 2 F-»A 2 F* . Again, there will be litde 
oscillation, so these generate a small component of shut 
times with a mean life close to that in AF. At low concen- 
trations, a channel that after opening reaches AR will 
very rarely rebind glycine and reopen, but this route 
will become increasingly common as glycine concentra- 
tion is increased, and it will contribute to the increase 
in duration of activations with concentration. 

Finally, the shut time between one activation and the 
next at a low concentration is predicted by the fit to be 
very long. At 20 uM glycine, the range is 25.6-92.9 s for 
the three sets used for fitting. The fact that it is so long 
explains why so few activations could be recorded, and 
it also explains why the P op en is so low (around 0.01) at 
20 uM glycine. The shut time between activations is 



characterized by at least one sojourn in the resting state, R. 
However, even at 20 uM glycine, the mean lifetime of 
a single sojourn in R is quite short, a few hundred milli- 
seconds. So, evidently, there is a lot of oscillation be- 
tween shut states between one activation and the next. 
Consider oscillations between the resting state, R (state 6) , 
and the monoliganded resting state, AR (state 5), and 
back. The distribution of the number of such oscil- 
lations can be found as follows. From R, the channel 
has to reach AR (state 5; i.e., tt 65 is 1). There will be no 
R-AR oscillations if the channel then moves straight to 
either state AF (state 3) or A 2 R (state 4) without returning 
to R. Therefore, the probability of there being no R-AR 
oscillations is P(0) = tt 53 + ir 54 = ( 1 — tt 56 ) • The probabil- 
ity of one oscillation (R->AR->R->AR->state 3 or 4) 
must be P( 1 ) = tt 56 ( 1 — tt 56 ) . So, in general, the proba- 
bility of r oscillations between R and AR before leaving 
for state 3 or state 4 is seen to be 

^) = (7t 5B ) r (l-7t 56 ). < 6 > 

These probabilities (from r= 1 to °°) sum to 1, and the 
mean number of oscillations is 

For example, in the fitted set that was used for 
the histograms, tt 56 = 0.996 at 20 uM glycine, so the 
mean number of R— >AR— >R oscillations is 229. The 
mean length of a sojourn in Rwas 251 ms, and the mean 
length of the monoliganded state AR is 2.17 ms. Thus, 
the mean length of a run of R— >AR— >R oscillations is 
229 (251 + 2.17) = 58 s. This is almost as long as the 
mean shut time between activations (74 s in this case) . 
Thus, the very long shut times that separate long activa- 
tions are due to the many R— >AR— >R oscillations rather 
than to the long life of any individual shut state. 
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